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We report a joint test of local Lorentz invariance and the Einstein Equivalence Principle for 
electrons, using long-term measurements of the transition frequency between two nearly degenerate 
states of atomic dysprosium. We present many-body calculations which demonstrate that the energy 
splitting of these states is particularly sensitive to violations of both special and general relativity. 
We limit Lorentz violation for electrons at the level of 10 -17 , matching or improving the best 
laboratory and astrophysical limits by up to a factor of 10, and gravitational redshift anomalies for 
electrons by two orders of magnitude, to 10 -8 . With some enhancements, our experiment may be 
sensitive to Lorentz violation at the level of 9 x 10 -20 . 
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Local Lorentz invariance (LLI) and the Einstein equiv- 
alence principle (EEP) are fundamental to both the stan- 
dard model and general relativity [TJ. Nevertheless, these 
symmetries may be violated at experimentally accessible 
energy scales due to spontaneous symmetry breaking, or 
some other mechanism at high energy scales [2] [3[ . This 
has motivated the development of many experimental 
tests of both LLI and EEP jH [5] , and of a phenomeno- 
logical framework, known as the standard model exten- 
sion (SME), which can be used to quantitatively com- 
pare these tests' results to one another [BJ. This widely 
used [1] framework augments the standard model La- 
grangian with every combination of standard model fields 
that are not term-by-term invariant under Lorentz trans- 
formations, while maintaining gauge invariance, energy- 
momentum conservation, and Lorentz invariance of the 
total action [6; . Violations of LLI, which themselves con- 
stitute violations of EEP [TJ [5] , have also been shown to 
violate other tenets of general relativity [7J. 

In this letter, we show, using many-body calculations, 
that the energies of two low-lying excited states of dys- 
prosium (Dy) [5141 lj are extremely sensitive to physics 
that breaks LLI and the EEP in the dynamics of elec- 
trons. We report the results of an analysis of Dy spec- 
troscopy data acquired over two years that significantly 
improves upon the best laboratory [12j and accelera- 
tor [13] limits on electron violations of LLI and EEP. Our 
result is competitive with some astrophysical bounds |14j . 
We also improve constraints on electron-related gravita- 
tional redshift anomalies |15j by two orders of magni- 
tude [16] . 

The EEP and LLI require that spacetime, while it may 
be curved, be locally flat, and Lorentzian pp. Thus the 
relative frequencies of any set of clocks at relative rest 
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and located at the same point in (or within a sufficiently 
small volume of) spacetime be independent of a) where 
that point is located in a gravitational potential, or b) 
the velocity of their rest frame (LLI). In the SME, vi- 
olation of EEP and LLI for electrons can be described 
by modifying the electron dispersion relation, which in 
turn causes the energies of bound electronic states to vary 
with the velocity and gravitational potential of their rest 
frame [7J H7]. 

We focus on the symmetric, traceless c M „ tensor in the 
electron sector of the SME, written using coordinates 
such that the speed of light is a constant c in all frames. 
The c^u tensor modifies the kinetic term in the electronic 
QED Lagrangian to become [BJ 

C = §ty (rf v + C/u ,7") D v ip— 4,m^j, (1) 

where m is the electron mass, ip is a four-component 
Dirac spinor, 7^ are the Dirac gamma matrixes, and 

fD v g= fD p g - gD v f , with D v = d v - qA v . The 
c^v tensor is frame-dependent [BJ I17H19) , and is uniquely 
specified by its value in a standard reference frame. We 
use the Sun-centered, celestial equatorial frame (SCCEF) 
for this purpose, indicated by the coordinate indexes 
(T, X, Y, Z), for ease of comparison with other re- 
sults [3]. The component indexes for laboratory frame 
coordinates are given as (0,1,2,3), where t = xq/c is 
the time coordinate. Roman indexes are used to in- 
dicate the spatial components of c^ v , and are capital- 
ized in the SCCEF frame. The Cp V tensor has six 
parity-even components: ctt, plus the five cjx's; and 
three parity-odd components: ctj, which introduce di- 
rection and frame dependent anisotropies in the elec- 
trons' energy-momentum, or dispersion relation [BJ. This 
shifts the energies of bound electronic states as a func- 
tion of the states' orientation and alignment in abso- 
lute space, breaking both LLI and rotational symme- 
try [T7J. In a gravitational potential, Eq. ([TJ acquires 
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FIG. 1. Energy levels of Dysprosium. Atoms are optically 
pumped (solid lines) to a state which decays (wavy lines) into 
the metastable state B. A linearly polarized rf field drives 
the B — > A transition, which is detected via fluorescence at 
564nm. Insets a) and b) show the magnified diagram for 
164 Dy and 162 Dy, respectively. Lorentz-symmetry violation 
shifts the rf resonance by 8lj t i = (AEb — AEa)/Ti. The sign 
of the observed shift depends on the sign of the level splitting. 

additional terms proportional to and to the curva- 
ture of spacetime [7J. These arise due to the interplay 
between the LLI-preserving distortion of spacetime due 
to gravity, and the LLI- violating distortion due to c p „, 
generating anomalous gravitational redshifts that scale 
with the electron's kinetic energy [71 ITS]. 

In terms of spherical-tensor operators, Eq. ([I]) pro- 
duces a shift Sh in the effective Hamiltonian for a bound 
electron with momentum p given by [7J 1171 120j 




where we have included the leading order (2[//3c 2 )coo 
gravitational redshift anomaly [7J [5T] in terms of the 
Newtonian potential U, and 

Cq ' 1 = c o + | c jii C^ 2 ' = ( c jj ~ 3c~33) 

Cfl = ±3(c 3 i ± ic 32 ), C { ±1 = 3(cn - c 22 ± icx 2 ), 

are written in terms of the laboratory-frame values of 
the c^u tensor, with summation implied over like indexes. 

Note that Cq is also known as c q in the literature [T7| . 
The spherical tensor components of the squared momen- 
tum are written as T^ 1 = p 2 — 3p|, = ±p 3 (p 1 ±ip2), 

and T±2 = {p 2 — P 2 )/2 ± ipiPi- The energy shift for a 
state I J, M) of an atom due to the perturbation |2| is 
the expectation value of the corresponding N electron 
operator. Since only tensors with q = contribute to en- 
ergy shifts of bound states, we need only calculate matrix 
elements for the p 2 and = p 2 — 3p 2 operators. 



TABLE I. Matrix elements of the relevant operators of 
Lorentz violation for the states A and B of Dy in units of 
the Hartree energy E h = (6.5 x 10 15 Hz)h. 



State A State B 

Term Symbol *\W] 7 H 5 

Energies (cm - ) 
Expt. [22] 19798 19798 

Calc. rjg| 19786 19770 

Matrix Element (units of Eh) 
(J || c 7 (7 j pj - 37 3 p 3 ) || J) 69.48 49.73 
(J || p 2 - 3pi || J) 69.84 49.89 

(JM\p 2 \JM) 437 422 



Dysprosium, an atom with 66 protons and a partially 
filled /-shell, is well suited to measuring the electron 
coefficients. It possesses two near-degenerate, low-lying 
excited states with significant momentum quadrupole 
moments, opposite parity, and leading configurations: 
[Xe]4/ 10 5d6s, J = 10 (state A) and [Xe]4/ 9 5ci 2 6s, J = 10 
(state B) 1 which differ by a transposition of an electron 
from the 4/ to the hd orbital. The energy difference be- 
tween these states can be measured directly by driving an 
electric-dipole transition (Fig. [IJ with a radio-frequency 
(rf ) field, and should be particularly sensitive to anoma- 
lies proportional to the electrons' kinetic energy, since the 
4/ orbital lies partly within the radius of filled s, p, and 
d shells that screen the nuclear charge from the larger 5d 
orbital. The same properties of Dy that make measure- 
ments of the B — > A transition particularly sensitive to 
variations in the fine structure constant, a [24Tl26j . also 
make them highly sensitive to violations of LLI and EEP. 

To calculate the relevant matrix elements for these 
states, we use a version of the configuration interac- 
tion method optimized for atoms with many electrons 
in open shells. This method has been used to calcu- 
late energy levels, transition amplitudes, dynamic po- 
larizabilities, "magic" frequencies in optical traps, and 
the effects of a variation and parity violation in Dy and 
other atoms [27j [28] • Calculated values of the reduced 
matrix elements for the A and B states of Dy are pre- 
sented in Table [TJ and details of their derivation can be 
found in the Supplement |29] . To check our results, we 
calculate the fully relativistic matrix elements of c~f°^pj 
and Tq 2 "* = C7°(7 J pj — ^l^Pz), corresponding to p 2 and 
p 2 — ip 2 [20]. We find good agreement between both 
calculations, consistent with our initial approximation. 

An effusive atomic beam of Dysprosium atoms is pro- 
duced by a ~ 1400 K oven, and is optically pumped into 
the metastable state B via consecutive laser excitations 
with 833 nm and 669 nm light, followed by a spontaneous 
decay. The atoms are resonantly excited from state B to 
A via an rf electric field, whose linear polarization de- 
fines the atoms' quantization axis. The polarization of 
the excitation laser is chosen to create a symmetric pop- 
ulation among the ±M magnetic sublevels of state B. 
The A state relaxes to the ground state in a cascade 
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FIG. 2. Full record of frequency measurements for 162 Dy 
(upper data set) and 164 Dy (lower data set). Frequencies are 
plotted relative to 234 661 065 Hz for 162 Dy and 753 513 695 
Hz for 164 Dy. Error bars are obtained by binning measure- 
ments into sets of 20 and calculating the standard error of the 
mean for each set. The solid line indicates the least-squares 
fit. Inset: an expanded view of the most recent measure- 
ments beginning on Oct. 19, 2012, with time given in Pacific 
Standard Time (Coordinated Universal Time minus 8 hours). 



TABLE II. Constraints on electron c^-coefficients from spec- 
troscopy of the rf transitions in 162 Dy and 164 Dy. We use 
the shorthand notation cx-y = cxx — cyy, c T (y+z) = 
cty cos 77 + ctz sin 77, and c T (y-z) = cty sin 77 — ctz cos 77, 
where 77 = 23.4° is the angle between the Earth's spin and or- 
bital axes. Bounds above the horizontal divider are obtained 
from 12 hours of continuous measurement, while those below 
the line are obtained from analysis of over two years of data, 
see text. Some uncertainties for the latter limits are adjusted 
for systematic error; the statistical uncertainty is then indi- 
cated in parenthesis. Past bounds on cjk, ctj, and ctt, 
and the purely gravitational limit on ctt are from analyses 
reported in [12], [14], [13], and [16], respectively. 



Combination 


New Limit 


Existing Limit 




0.10 cx-y + 0.99 c X z 


-9.0 + 11 


-33 + 19 xlO" 


17 


0.99 cx-y - 0.10 cx z 


3.7 + 5.6 


-8.6 + 62 xl0~ 


17 


0.94 c X y + 0.35 c Y z 


-0.4 + 2.8 


43 + 19 xl0~ 


17 


0.35 c X y - 0.94 cyz 


3.2 + 7.0 


5.3 + 23 xl0~ 


17 


0.98 c T (y+z) — 0.19 ctx 


-1.3 + 18(3.3) 


0.3+1.8 xl0~ 


15 


0.19 c T (y+z) + 0.98 Ctx 


7.3 + 6.6(2.4) 


-1.5 + 5.4 xl0~ 


15 


CT(Y-Z) 


-38 + 30(5) 


.011 + .028 XlO" 


13 


ctt 


-2.1 ± 6.7 


10~ 6 (2±2) xl0~ 


9 


Ctt (gravitational) 


-14 + 28(9) 4600 + 4600 xl0~ 


9 



decay, emitting 564 nm light in the process. The transi- 
tion frequency is determined by measuring the intensity 
of the 564 nm fluorescence (with a photomultiplier) as a 
function of radio frequency. More details regarding the 
experimental procedure and apparatus can be found in 
Refs. [25] [30]. 

We measure the average frequency shift of all popu- 
lated magnetic sublevels of state B relative to those of 
state A that are coupled by the rf electric field. The en- 
ergy shift of each transition is calculated using Eq. (H 
and the calculated reduced matrix elements for each 
state. The actual distribution of population among the 
magnetic sublevels is found by resolving the Zeeman 
structure of the two states, and measuring the peak am- 
plitude of each transition. These amplitudes are used 
as weights in a sum of the shifts of each state due to 
Eq. ([2]) to determine the average shift of the unresolved 
line. The average shift in the B — > A transition frequency 
ujtf is given by 



^ = ±(10 14 Hz) 



r>ou|cf-^coo 



9.1C, 



(2: 



, (3) 



where Uq = —MQG/ri^b is the Sun's gravitational po- 
tential, and Wrf is defined to be positive, producing a 
positive (negative) shift for 164 Dy ( 162 Dy). This sign dif- 
ference helps reject background systematics, and is de- 
termined by that of the splitting between A and B. The 
sign of the second term depends on the relative magnetic 
sublevel populations. 



The value of Cq ^ and C^f' in the laboratory-frame 
is a function of Cau in the SCCEF, and the orientation 
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and velocity of the lab. Thus any anomalous duj r f mea- 
sured in the lab must vary in time [171 . The precise 

relation between and Cg 2 ' and the SCCEF value of 
can be found in the Supplement [55] . The scalar Ctt 
component of can be bounded via frame- or gravi- 
tational potential-dependent effects, as it contributes to 

(2) 

the modulation of Cq , scaled by the Earth's orbital ve- 
locity squared ^»lx 10~ 8 , and to that of the larger 
scalar term in Eq. ^ via modulations of the laboratory in 
the Sun's gravitational potential, which have amplitude 

au q /c 2 = 1.7 x nr 10 . 

Using repeated measurements of Slu^ acquired over 
nearly two years, we obtain constraints on eight of the 
nine elements of c M „.The cjk coefficients are constrained 
using data collected over the course of 12 hours beginning 
on Oct. 19, 2012. For each isotope the mean value of 20 
successive frequency measurements (~ 10 seconds) is as- 
signed an error bar according to the standard error of the 
mean for that bin. The resulting data are fit to Eq. ([3| in 
terms of cjk in the SCCEF [29] . augmented by an inde- 
pendent, constant frequency offset for each isotope. The 
short duration of this data set allows us to neglect the 
slow (1 yr _1 and 2 yr _1 ) variations induced by the ctt 
and ctj terms. These terms are neglected in this fit, as 
they are suppressed by at least one factor of /3 ffi ~ 10~ 4 , 
and existing limits P3j on these terms constrains their 
contributions well below our statistical sensitivity. 

The ctj and ctt coefficients are constrained using 
data collected between November 2010 and July 2012. 
The data are binned and assigned error bars as previ- 
ously described. Since the above analysis of the 12 hour 
data set provides tight constraints on cjk coefficients, 



4 



5.0 
^ 2.0 

N 

£ i.o 

£ 0.5 

b 

0.2 
0.1 



' 1 C33 ^ 2^ X 1 _ _ _ '••».♦... if 



1 



5 10 50 100 

Averaging time, t (s) 



, Is 

500 1000 



FIG. 3. Allan deviation from a two hour measurement of the 
164 Dy transition frequency (7:00 to 9:00 on the Fig. [2] inset). 



the second fit includes only the ctj and ctt coefficients. 
The fit routine is the same as before, adding an indepen- 
dent linear slope for each isotope to account for long-term 
systematic drifts. The resulting fit includes a large signal 
for the combination ct(y-z) = c ty sin rj — ctz cos ?y = 
(—38 ± 5) x 10~ 13 , where 77 is the Earth's axial tilt. 
As such a signal is inconsistent with existing limits on 
Ct(y-z) [SI EH: we suspect the presence of uncontrolled 
systematic shifts in 5u) r { with characteristic modulation 
frequencies near 1 day -1 and 2 day -1 , and amplitude 
300 mHz. These systematics may be due in part to, e.g. , 
magnetic field fluctuations (~ 50 mHz), blackbody shifts 
due to changes in the temperature of the spectroscopy 
chamber (~ 60 mHz) [31] , and changes in electronic off- 
sets (~ 140 mHz). Daily fluctuations in these system- 
atic shifts have less effect on our bounds on ctx and 
c t{y+z) = cty cos 77 + ctz sin 77, as these are primarily 
sensitive to the yearly modulation signal produced by the 
larger scalar component of Eq. ^ [29] . In the presence 
of correlated noise, our statistical error bars overstate 
our measurement's precision, thus we repeat the least- 
squares analysis without flipping signs for 162 Dy rela- 
tive to 164 Dy. This model is insensitive to Lorentz viola- 
tion, but is sensitive to systematic error. Where they are 
larger, the absolute mean of each term in this fit replaces 
the statistical error estimated by the original fit. 

We have also analyzed our results as a test of the grav- 
itational redshift for electrons in the Sun's gravitational 
potential by fitting the long term data to terms propor- 
tional to the gravitational potential, neglecting frame- 
dependent effects. We obtain a purely gravitational limit 
on the electron's ctt coefficient of —14 ± 28 x 10 -9 . 

The data and fits are shown in Fig. [2j The fit results 
are displayed in Table [H] with uncertainties quoted for 
68% confidence limits. The reduced chi-squared, x 2 , for 
the short and long time-scale fits are 1.2 and 1.8, respec- 
tively. The larger % 2 of the long-term fit is likely due to 
uncontrolled systematics that have not been accounted 
for in our purely statistical estimation of error bars. To 
obtain conservative estimates on parameter uncertainties 



we have scaled the statistical error bars in both fits to 
provide \ 2 = 1. 

We have tightened experimental limits on four of the 
six parity-even components of the c M „ tensor by factors 
ranging from 2 to 10 |U [12] ■ We report limits on two 
combinations of the parity-odd ctj that are on par with 
those set by ~ 50 TeV astrophysical phenomena [H Q3] . 
We improve bounds on electron-related anomalies in the 
gravitational redshift by a factor of 160, to 2.8 x 10~ 8 . 
With optimization, our experiment could yield signifi- 
cantly improved constraints. As Fig. [3] shows, our ex- 
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periment is statistically sensitive to Cq 



3c. 
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in the lab at the level of 2.2 x 10~ 16 after 400 seconds 
of averaging. At present, we must wait a full day for 
the Earth to rotate the laboratory in the fixed reference 
frame, increasing our susceptibility to systematics vary- 
ing on that timescale. This could be addressed by active 
rotation of the entire apparatus, or of the polarization 
of the rf electric field, making possible statistically lim- 
ited sensitivities to at the level of 1.5 x 10~ 17 in one 
day, and 7.8 x 10 -19 in a year. Optically pumping the 
atoms to the M — ±10 states could increase the experi- 

(2) 

ment's sensitivity to Q by a factor of ~ 4.5. Increas- 
ing the interaction time of the atoms in the rf field could 
gain another factor of two, as the measured linewidth 
of 40 kHz is twice the natural linewidth of state A. An 
optimized experiment may thus reach sensitivities at the 
order of 8.7 x 10 -20 in one year. This would be three 
orders of magnitude better than the presently reported 
limits on cjk, two orders of magnitude better than the 
best sensitivities attainable by existing optical resonator 
tests |32j , and could prove more sensitive than astrophys- 
ical tests O [33] . Still narrower linewidths are possible 
in spectroscopic measurements of the Zeeman and hy- 
pcrfine structure of the ground state of trapped Dy [54"] . 
other rare-earth elements, and of the long-lived states of 
rare-earth ions in doped materials. Optical transition en- 
ergies in trapped ion or neutral atom clocks, and of the 
electronic and ro- vibrational states of molecules may also 
be sensitive to c M „. The latter might also probe violations 
of LLI and EEP for nuclei. The derivation of the scalar 
and quadrupole moments of the involved states will be 
the subject of future work. 
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SUPPLEMENTARY MATERIALS: 
ELEMENTS 



MATRIX 



Here, we present the derivation of the explicit form 
of the single orbital matrix elements for the operators 



p 2 and T, 



(2) 



— y — 3^3. We also present the matrix 

-(2) _ 



elements for the relativistic operators c-f°^Pj and T 

c 7°(7 J Pj ~ 37 3 p 3 )- 

The many-electron state of an atom is a linear com- 
bination of Slater determinants, constructed from single- 
electron orbitals. These orbitals, in the relativistic limit, 
can be written as 



1 



r \iag nK (r)Cl. 



(4) 



where n = =f(j + 1/2) (with j = I ± 1/2) is the quantum 
number for the angular momentum of a Dirac spinor, 
a is the fine structure constant, and the spin-dependent 
component Q Km (9,<p) is given by 



ft* 



± K+ i/2- m y 



2k+1 



r l,m-l/2\"i 1 



K+l/2+m 
2k+1 



Y, 



(5) 



where Y^ m (9, <$>) is a spherical harmonic. The upper com- 
ponent of Eq. Q corresponds to the electron component 
of the spinor wave function, while the lower component 
to the positron. In the non- relativistic limit, we can drop 
the antimatter component of Eq. Q and evaluate matrix 
elements in terms of 



f-* Inn \1j q 



(6) 



Using the Pauli matrix identity (<? • p) 2 = p 2 , with a a 
vector of Pauli matrixes, we can write the action of the 
operator p 2 on ip nK , m (r) as 



{S ■ p) 2 tp nKm {r) = 



nr f d 2 f nK _ k(k + i) 

dr 2 r 2 



friK, I 



(7) 



() 



The matrix elements for p 2 are thus equal to S K ' K Ii, 
where I± is the radial integral 



h = hf / dr 
Jq 



dfn'K- df nK k(k + 1) 
dr dr r 2 



fn'tz' fn 



, (8) 

Similarly, the matrix elements of p\ = —Ti d 2 /dz 2 can 
be found by applying p^ = —ifid/dz to the spinor twice 
{e.g., by using Eq. (3.3.3) of [35]). Using the Wigner- 

(2) 

Eckart theorem, the matrix elements of T can be writ- 
ten in terms of reduced matrix elements after factoring 
out the 3j-symbol: 



f 2 j 

-TO TO 



(nVTO|T (2) |wc7n) = {-iy'- m 

x (nV||T (2) ||nK). (9) 
These reduced matrix elements are then given by 



(n'K'\\T^\\n K ) = < 



\A(f,j)h ifn 

A(j',j)h if« 

B(j'J)I 2 if ^ 

B(j',j)I 3 ifn 

[C(j)h if« 



= — K — 1 

= -K + 1 

= K — 2 

= K + 2 

= K 



(10) 



where A(j',j), B(j'J), and C(j) are given by 



6(2j m + 3)(2j m + l) 
(2j m + 2)2j m (2j m +4) 



1/2 



B(f, j), and C(j) are the same, while the radial integrals 
differ, so that 

(n' K '|| C7 ( 7 ^- 37 3 p 3 )|^) 

( At Al A\ ' 



if k' 


= — K — 


if k' 


= -K + 


if k' 


= K-2 


if k' 


= K + 2 


if k! 


= ft 



where 

A = - 



dr ( (2k - l)g n 'K'^^ 
o V dr 



(2 K + 3)f n , K . 



car 
"~2~ 

(2k -!)(« + !) (2k + 3)k 

r r 



dr 



car 
~~2~ 



drl (2k 3://,. , • (2k ■ 1 :/,..,.• 



dr 



dr 



(2«-3)* „ (2k + !)(«-!) , 



h = 2cah I drl gn'^'H^ 1 — Snv/nA 

I 5 = caTi j dr I g„/ K / — / n / K / 



r r 



3(2j ro + 5)(2j m + 3)(2j m + l) 
2(2j m +4)(2 Jm + 2) 



1/2 



SUPPLEMENTARY MATERIALS: MANY-BODY 
CALCULATION 



C(i) = 



(2j + 3)(2j + l)(2j-l) 
4j(j + 1) 



1/2 



where j m — min(j', j). The radial integrals I 2 and ^3 are 



I 2 =fi dr. 

/□ V or or 



9/n' K ' df nK 2k - 1 d/„ K 



r " Or 
k(k - 2) 



fn'K' fnn 



and 



/o V <'>' '''' 



9/n' K ' <9/nK , 2ft + 3 d/„ K 



/ n' k 

r Or 



(k + 3)(« + 1) 



Jti'k,' fnK 



To evaluate the expectation of the electronic kinetic 
energy for states A and B in Dy, we use a version of the 
configuration interaction (CI) method, originally devel- 
oped for calculating energy levels, electromagnetic am- 
plitudes, dynamics polarizabilities, "magic" frequencies 
in optical traps, the effects of variations in the fine struc- 
ture constant, and parity violation in atoms with many 
electrons in open shells [37]. The effective Hamiltonian 
for N v valence electrons in an atom (N v = 12 for Dy) 
has the form 



H 



cff 



N v 



E< 



( r i) + J2 e2 / r ik, 



(12) 



where e is the elementary electron charge, and h\j is the 
one-electron Hamiltonian 

Ze 2 

h Xj {r 5 ) = ca-p+(fi- l)m e c 2 + V core + SV, (13) 



As for the the matrix elements of the relativistic form 
of the operators acting on the four-component spinor in 
Eq. Q, we find that their angular components A(j',j), 



where a, = 7°7j, and j3 — 7 . Here, V core is the Hartree- 
Fock potential due to the core electrons, and SV simu- 
lates the effects of correlations between core and valence 
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TABLE III. Configurations and effective core polarizabilities 
a p (a.u.) used in the calculations. 



N 


Parity 


Configuration 


0!p 


1 


Even 


4/ lu 6s 2 


0.4 


2 


Even 


4/ 10 6s5d 


0.4006 


3 


Even 


4/ 9 6s 2 6p 


0.4039 


4 


Even 


4/ 9 5d6s6p 


0.389 


5 


Even 


4/ 10 6p 2 


0.4 


6 


Even 


4/ 9 5d 2 6p 


0.4 


7 


Odd 


4/ 9 5d 2 6s 


0.3947 


8 


Odd 


4/ 9 5d6s 2 


0.3994 


9 


Odd 


4/ 10 6s6p 


0.397 


10 


Odd 


4/ 10 5d6p 


0.4 


11 


Odd 


4/ 9 5d6p 2 


0.4 


12 


Odd 


4/ 9 6s6p 2 


0.4 



electrons. It is also known as the polarization potential, 
and has the form 



SV = 



2{r 4 



(14) 



where a p is the core polarizability, and a is a cutoff pa- 
rameter (here, the Bohr radius a = as)- Table III lists 

The 



the configurations considered in our calculation, 
self-consistent Hartree-Fock procedure is performed sep- 
arately for each configuration. Next, valence states ob- 
tained from the Hartree-Fock calculations are used as ba- 
sis states for the CI calculation. The CI method requires 
that the atomic core remain the same for each outer elec- 
tron configuration. We select the core state that corre- 
sponds to the ground state configuration. Changes in the 
core state due to changes in the valence state are small, 
and can be neglected, as the 6s, 6p and 5d shells are 
comparatively distant, and the potential they produce at 
the core is a nearly uniform perturbation. We can thus 
model the valence potential as shifting the energy of the 
core eigenstates, without modifying the wave functions. 
The 4/ electrons, on the other hand, lie nearer to the 
core, and thus have a stronger effect. In all cases con- 
sidered here, however, only one of the ten 4/ electrons 



ever change state. Thus their total effect on the atomic 
core will also be small. A more detailed discussion of the 
effects of valence electrons on atomic core states can be 
found in Refs. M. 



The form of SV in Eq. ( 14 ) is chosen to coincide with 



the standard polarization potential at large distances 
(— a p /2r 4 ). We treat the a p for each configuration as a fit 
parameter, chosen to match experimentally measured en- 
ergy intervals between states of different configurations. 
Their values are displayed in Table |III| The values of 
a p are similar for all configurations considered. This is 
expected, as the state of the core electrons is the same 
for each valence configuration. The small differences in 
a p for different configurations help compensate the effect 
of the incompleteness of our valence electron basis, and 
other imperfections in the calculation. 

SUPPLEMENTARY MATERIALS: FRAME 
DEPENDENCE OF Suj ri 

The laboratory frame components of the c M „ tensor are 
found in terms of the sun-centered equatorial frame (SC- 
CEF) components by requiring that the Lagrangian of 
Eq. (1) be invariant under the observer Lorentz transfor- 
mation between frames. This condition leads to 



A M A N 



(15) 



where cmn is the c tensor defined in the SCCEF. In the 
SCCEF the Z axis points along the Earth's rotation axis, 
the X — Y plane is the Earth's equatorial plane, and the 
Y axis points from the Sun to the Earth at the vernal 
equinox. The Lorentz transformation A^ is a rotation, 
dependent on experimental geometry and the rotation of 
the Earth, to align the laboratory defined axes with these 
SCCEF axes followed by a boost, determined mainly by 
the Earth's orbital velocity, to the rest frame of the Sun. 
Our laboratory frame is defined such that z is parallel to 
the quantization axis defined by the polarization of the 
rf field. Thus z points 9 — 15° south of west, x points 
9 = 15° east of south, and y points vertically upwards. 
The rotation from the SCCEF to the laboratory frame is 
given by 



cos x cos 9 cos w7® — sin 9 sin wT® sin 9 cos wT© + cos x cos 9 sin wT© — sin x cos 6 N 
R = | — sin x cos ujT® — sin x sin wT® — cos x 

^— cos x sin # cos wT© — cos (9 sin cjT© cos 9 cos wT© — cos x sin 9 sin wT© sin x sin t 



(16) 



where x = 52.1° is the colatitude of our laboratory, ui 
is the angular frequency of the Earth's rotation in the 
SCCEF (i.e. 2irx an inverse sidereal day), and T® is 
measured from the first time that East, as measured in 
the laboratory, and the SCCEF Y-axis coincides after 
a vernal equinox. The boost of the laboratory in the 



SCCEF frame is given by [4] 

(/3®sinf2T — (3 l sin x sin. ojT \ 
- (3® cos ?] cos £IT + f3 L sin x cos ojT , (17) 
— /3® sin rj cos flT J 

where rj = 23.4° is the angle between the eclip- 
tic plane and the Earth's equatorial plane, /3© = 



TABLE IV. Dominant time-varying terms in the fit for the cjk coefficients. The frequencies a;® and O and the sidereal-day and 
yearly frequencies, respectively. The colatitude of the experiment is given by \ ~ 52.1°, 9 ~ 15° is the angle the quantization 
axis is rotated towards the South from West, and 77 ~ 23.4° is the angle between the ecliptic and the Earth's equatorial plane. 
The orbital boost is /3® ~ 10" 4 . The constants S ~ +5 x 10 16 Hz and Q ~ +9.1 x 10 14 Hz are the scalar and quadrupole shifts, 
respectively, from Eq. pi. For ease of comparison with Table II, we have defined cx-y = cxx — cyy and cx+y = cxx + cyy ■ 
Additional terms of 0{Qj3%) have been suppressed, although they are included in our fits. 









w® 

2w® 

20 


3Q (cxz sin x sin 28 + cyz sin 2\ sin 2 9) 
-|S [cx-Y cos x sin 29- \c X y (l + 3cos29 - 2cos2xsin 2 9)] 


3Q (cx z sin 2x sin 2 9 — cyz sin x sin 29) 

r cxr cos x sin 29+ 
^ ^cjc-y (l + 3 cos 29 -2 cos 2 X sin 2 9) 

5 £-» a2 • 2 



TABLE V. Dominant time-varying terms in the fit for the ctj and ctt coefficients. The frequencies w® and O and the 
sidereal-day and yearly frequencies, respectively. The colatitude of the experiment is given by x ~ 52.1°, 9 ~ 15° is the 
angle the quantization axis is rotated towards the South from West, and 77 ~ 23.4° is the angle between the ecliptic and 
the Earth's equatorial plane. The orbital boost is /3® ~ 10~ 4 . The constants S ~ +5 x 10 16 Hz and Q ~ +9.1 x 10 14 Hz 
are the scalar and quadrupole shifts, respectively, from Eq. pi. For ease of comparison with Table II, we have defined 
Ct(y+z) = cty cos 77 + ctz sin 77 and c T (y~z) = cty sin 77 — ctz cos 77. Additional terms of 0(<2/3®) proportional to ctx and 
Ct{y+z) have been suppressed, although they are included in our fits. The gravitational terms appear at sinfiT and cosQT due 
to a phase offset <j>Q = 10.4° between the oscillation of the boost vector (measured from the vernal equinox) and the oscillation 
of the Earth in the Solar gravitational potential (with perihelion on Jan. 3rd). The amplitude of the Earth's modulation in 
the Solar gravitational potential is AU/c 2 w 1.7 x 10 -10 . 

CO 7 Si C-i 




20 

w® — 20 
w® — O 
w® 

w® + O 
uj® + 20 

2wg 
2u7g 
2u7g 

2oj s 

2uJa 



- 20. 

- o 

+ Q 
+ 20. 



fSP 9 CTX + ^S-Sctt sin 



2QP%c TT cos 3 § sin § sin 2 9 sin 2x 
|Q/3ffiC T{ y_ z) (cos 77 + cos2?7) sin 2x sin 2 9 

- \ Q/3® ctt sin 2r) sin 2x sin 2 9 

— 1 2/?®<>r(-r-z) (cos 77 — cos 2rj) sin 2x sin 2 « 
\ QP%ctt (2 sin 77 — sin 2?7) sin 2x sin 2 9 
QP% ctt cos 4 I cos x sin 29 
|Q/?ffiC T (y_z) sin 77(1 + cos 77) cos x sin 29 

— \ QP%ctt sin 2 77 cos x sin 29 
3Q/?ffiC T (y_z) sin 77 sin 2 | cos x sin 29 
Q/3®ctt sin 4 2 C os x sin 29 



-fS/3 e CT (Y+ Z) + ^Sctt cos 0Q 
+ |Q/3®CT(y-z) sin2?7 (l + 3 cos 29 + 6 cos 2x sin 2 ( 
§ QP%c TT sin 2 77 (1 + 3 cos 29 + 6 cos 2x sin 2 9) 

— \ QP% ctt sin 77 ( 1 + cos 77) sin x sin 29 

— 2 Q/ 3 ® c T (y _ z) (cos 77 + cos 277) sin x sin 29 

— \ QP% ctt sin 277 sin \ sin 29 
1 Q/3ec T (y_z) (cos 77 — cos 277) sin x sin 29 



sin 277) sin x sin 29 



\QPIctt (2 sin 77 

- 1 QP%ctt cos 4 2 (1 + 3 cos 29 - 2 cos 2x sin ; 
1 QP®Cr(Y-Z) 
g 2/3® ctt sin 77 
-|Q/3®CT(y- 



sin 77 (1 + cos 77) (l + 3 cos 29 — 2 cos 2x sin 2 9) 
' 1 + 3 cos 29 - 2 cos 2x sin 2 9) 
in 77 sin 2 § (l + 3 cos 29 - 2 cos 2x sin 2 9) 
. , , sin - ^ (l + 3 cos 29 - 2 cos 2x sin 2 9) 



2n(l a.u.)/c(l yr) ~ 10~ 4 is the boost from the Earth's 
orbital velocity, /5l = i?ew ~ 1.5 X 10~ 6 is the boost from 
the rotational velocity of a point on the Earth's equator, 
27r/f2 is a sidereal year, and T is the time since the epoch, 
chosen to be the vernal equinox in the year 2000 [I] . The 
boost /3r,sinx, though included in our fits, is too small 
to make a significant contribution to our fits, and will 
henceforth be dropped from this discussion. 

Since the transformation between frames is time- 
dependent, constant values of cmn in the SCCEF give 
rise to time varying frequency shifts in the laboratory 
value of c^ u , and thus to time variations in &j r f, via 
Eq. Since cmn also gives rise to an anomalous 

gravitational redshift, there is an additional contribu- 
tion proportional to ctt, such that 6u r { ps + (5 x 10 16 
Hz) x 2AC//(3c 2 ) cos(flT— ), as given in Eq. ([3]), where 
AU ~ 1.7 x 10 _10 c 2 is the amplitude of the Earth's yearly 



modulation in the solar gravitational potential due to 
the eccentricity of the Earth's orbit, and 4>& is such that 
cos(£!T — 00 ) is minimized at perihelion (^Jan 3). 

The cjk coefficients contribute leading order energy 
shifts at much shorter time scales (daily) than the ctj 
and ctt coefficients (yearly). As such, the cjk coeffi- 
cients are constrained using a single long data set ac- 
quired over the course of one day to minimize the influ- 
ence of systematic effects acting on long time scales. To 
fit the Cjk terms, we assume that the parity-odd Ctj 
coefficients are as constrained by astrophysical observa- 
tions [14] , and that the contribution of the ctt term to 
the daily modulation of 8u] T { is negligible. Although we 
retain all terms up to 0((3q) in our fit function, the dom- 
inant terms that are relevant to our fit are 

6oj Ti = A + ^ (Sj sin UjT + C 3 cos cjjT) , (18) 

i 
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where A is an independent offset parameter for each iso- 
tope, and the relevant frequencies coj, and quadratures Sj 
and Cj which contain the relevant components of Cjk, 
are summarized in Table |IV| Since our experiment al- 
ternates between measuring 5uj t [ for 162 Dy and 164 Dy, 
we do not directly compare 5u} T f for the two isotopes. 
We therefore perform a simultaneous fit of the separate 



functions ( 18 ) for each isotope, subject to the constraint 



that the values of cjk must be the same for both. 

The ctj and ctt coefficients are constrained using 
data acquired at irregular intervals over two years. To 
fit the ctj and ctt terms, we set cjk to zero, and drop 
terms proportional to /3m that do not multiply ctt, since 
existing bounds on cx+y = cxx + cyy are sufficient to 
ensure its doubly boost-suppressed contribution to duirf 
is negligible at our level of precision [3] . The fit function 
is 

Suiri = A + MT + {Sj sin uj,T + Cj cos uijT) , (19) 



where as before, A is an isotope-dependent offset, and the 
MT term is applied to remove any linear drifts. The rel- 
evant frequencies u>j, and quadratures Sj and Cj, which 
contain the relevant components of ctj and Ctt, are sum- 
marized in Table |v} We perform a joint fit on the 162 Dy 
and 164 Dy data as before, and note an apparently large 
signal for ct(y-z) ~ (38 ± 5) x 10 -13 , roughly 7.6 times 
the statistical error bar. As such a result would be in- 
consistent with other experiments [1] , we suspect it may 
be due to modulated systematic errors that we cannot 
fully distinguish from our model function with the exist- 
ing dataset. To estimate these errors, we fit the same 
data to a modified fit function that does not change sign 



for the two Dy isotopes, thus maximizing our sensitivity 
to any common-mode systematics, and replace our orig- 
inal statistical error bars on each coefficient with their 
magnitudes in the fit to the modified function, if they 
are greater. If our original signal for ct(y-z) were not 
due to systematics, and thus truly due to Lorentz sym- 
metry violation, we would expect to obtain a small value 
for Or(Y-z) m the fit to the modified function. Instead, 
we observe that the mean value of ct(y-z) without iso- 
topic sign reversal is nearly as large as it is for the fit to 
the Lorentz- violating model, and thus conclude that this 
is not evidence for violation of LLI. On the other hand, 
our estimates of the systematic error in our fits to ctx 
and ct(y+z) is comparatively quite low. This is consis- 
tent with the hypothesis that our systematic background 
modulates with a period of one or half a day, and aver- 
ages out over longer times, since our sensitivity to ctx 
and ct(y+z) comes primarily from the once-yearly mod- 
ulated scalar term in Eq. (J3j|, while Ct(y-z) is equally 
sensitive to both yearly and daily modulations propor- 
tional to the quadrupole term (see Table |V| . 

In both cases, the uncorrelated combinations of coef- 
ficients reported in Table II were found by diagonaliz- 
ing the covariance matrix from the least-squares fit to 
Eqs. (pi) and (fl9|). 



All frequency measurements are made relative to a 
Cs frequency reference. Though in general, violation of 
LLI may affect the frequency of this reference, such ef- 
fects can be neglected here. The relevant quadrupole 
term is strongly suppressed due to the symmetries of 
the 6s 1 state [17], and the scalar shift is expected to be 
~ (10 GHzjc,,,; [IB], while as Table I shows, the splitting 
between the A and B states of Dy will be ~ (100 THz)c M „ 
or greater. 



